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Abstract. Refocusing, or dynamical decoupling, is a coherent control technique where the internal dynamics of 
a quantum system is effectively averaged out by an application of specially designed driving fields. The method 
has originated in nuclear magnetic resonance, but it was independently discovered in atomic physics as a 
"coherent destruction of tunneling" . Present work deals with the analysis of the performance of "soft" refocusing 
pulses and pulse sequences in protecting the coherence of a qubit system coupled to a quantum oscillator. 



Introduction and background. Quantum coherent con- 
C^l trol has found way into many applications, including nuclear 
£H magnetic resonance (NMR), quantum information process- 
^ing (QIP), spintronics, atomic physics, etc. The simplest 
' 4 Control technique is dynamical decoupling (DD), also known 
as refocusing, or coherent destruction of tunneling. The 
g oal here is to preserve coherence by averaging out the un- 
c~| wanted couplings. This is achieved most readily by running 
O precisely designed sequences of specially designed pulses [1] . 
+1> , In a closed system this can be analyzed in terms of the 
Ch average Hamiltonian theory in the "rotating frame" defined 
r^by the controlling fields [2]. To leading order, the net evo- 
^4ution over the refocusing period is indeed described by the 
'—Hamiltonian of the system averaged over the controlled dy- 
| hamics. For an open system, the dynamics associated with 
^the bath degrees of freedom can be also averaged out, as 
llbng as they are sufficiently slow. This can be understood 
^^by noticing that the driven evolution with period r = 2tt/H 
*0 shifts some of the system's spectral weight by the Floquet 



-, harmonics, uj 



uj - 



nfi. With the average Hamiltonian for 



the closed system vanishing, the original spectral weight at 



n = disappears altogether, and the direct transitions with 
the bath degrees of freedom are also suppressed as long as 
• f2 exceeds the bath cut-off frequency, fl ^ oj c [3] . 
, ^ For a closed system, the refocusing can be made more ac- 
curate by designing higher-order sequences, where not only 
5— (the average Hamiltonian (k = 1), but all the terms of order 
^ k < K in the Magnus (cumulant) expansion of the evolution 
operator over the period r are suppressed. The correspond- 
ing simulation can be done efficiently by constructing time- 
dependent perturbation theory on small clusters [4]. The 
quantum kinetics of the corresponding open system with 
order- if DD, K < 2, was analyzed by one of the authors 
using the non-Markovian master equation in the rotating 
frame defined by the refocusing fields [5]. This involved a 
resummation of the series for the Laplace-transformed re- 
solvent of the master equation near each Floquet harmonic, 
with subsequent summation of all harmonics. 

The results of Ref. [5] can be summarized as follows. 
With K > 1 refocusing, there are no direct transitions, 
which allows an additional expansion in powers of the small 
adiabaticity parameter, ui c /Q. In this situation the decoher- 



ence is dominated by reactive processes (dephasing, or phase 
diffusion). With K = 1, the bath correlators are modulated 
at frequency f2. This reduces the effective bath correlation 
time, and the phase diffusion rate is suppressed by a fac- 
tor cx oj c /CI. With K = 2 refocusing, all 2nd-order terms 
involving instantaneous correlators of the bath coupling are 
cancelled. Generically, this leads to a suppression of the 
dephasing rate by an additional factor oc (w c /f2) 2 , while in 
some cases (including single-qubit refocusing) all terms of 
the expansion in powers of the small adiabaticity parameter 
(u) c /Q) disappear. This causes an exponential suppression of 
the dephasing rate, so that an excellent refocusing accuracy 
can be achieved with relatively slow refocusing, Q^j0j c . 

The conclusions in Ref. [5] were based on the analysis of 
the oscillator bath with a featureless spectral function, with 
the cut-off frequency oj c serving as the only scale describing 
the bath correlations. They do not apply in the presence of 
sharp spectral features which appear if the controlled sys- 
tem is coupled to a local high-Q oscillator. On the other 
hand, the situation where the controlled system is coupled 
to a local oscillator mode is quite common. This situation is 
realized in atomic physics, where the oscillator in question 
is the cavity mode, while the continuous-wave (CW) exci- 
tation is used to suppress the coupling. In several quantum 
computer designs, nearly-linear oscillator modes are inher- 
ently present (e.g., mutual displacement in ion traps, or QCs 
based on electrons on helium) . Finally, there are suggestions 
to include local high-Q oscillators in the QC designs to serve 
as "quantum memory" or "quantum information bath" [6]. 

In this work we consider refocusing of a qubit system 
where the spectral function of the oscillator bath has a sharp 
resonance. More specifically, we include the resonant mode 
in the system Hamiltonian, and consider the quantum kinet- 
ics of the resulting system in the presence of a featureless 
low-frequency oscillator bath driven by the refocusing pulses 
applied to the qubits only. Such a system can be analyzed 
with the help of the general results [5], as long as one is 
able to construct a K = 1 or K = 2 refocusing sequence 
to decouple the oscillator and other degrees of freedom. To 
this end, and having in mind sequences of soft pulses, we 
consider the analytical structure of the evolution operator 
for a closed system of arbitrary complexity, where one of 
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the qubits is driven by a single one-dimensional 7r-pulse. 
An analysis of any refocusing sequence is then reduced to 
computing an ordered product of evolution operators for in- 
dividual pulses. We illustrate the technique by analyzing 
the controlled dynamics of a qubit coupled to an oscilla- 
tor. One of the constructed sequences provide order K = 2 
qubit refocusing for any form of qubit-oscillator coupling, 
and was also shown to provide an excellent decoupling in 
the presence of a thermal bath. 

Single 7r-pulse. Consider a qubit with generic couplings, 

H s = <t x A x + a y A y + a z A z + A , (1) 

where cr M are the qubit Pauli matrices and A v , v = 0, x, y, z 
are the operators describing the degrees of freedom of the 
rest of the system which commute with cr^, [oy,,^] = 0. 
The qubit evolution is driven by a one-dimensional pulse, 

H c = & x V x {t), 0<t<T p , (2) 

where the field V x (t) defines the pulse shape. The evolution 
due to the pulse is dominant; the unitary evolution operator 
to zeroth order in Hs is simply 

U (t) = e"^*)/ 2 , </>(t) = f dt' V x (t'). (3) 

Jo 

When acting on the spin operators, this is just a rotation, 
e.g., Uo{t)<j y Ul{t) = a y cos (f)(t) + a z sin (j){t). Suppose V x (t) 
be symmetric, V x (t p — t) — V x (t), 7r-pulse, 4>{t p ) — ir, and it 
additionally satisfies the first-order self-refocusing condition 
s = (smcp(t))p = 0, where (f(t)) p denotes the time-average 
over pulse duration. Then, the evolution operator X = 
Uo(t p ) expanded to second order in t p Hs reads 

X ( V = -i ax - r p (A x + a x A ) + ^t 2 {A q , A x } 

+l^a x (A 2 + A 2 X ) + T 2 a(A 2 y + A 2 Z + ia x [A y , A z ]) 

+T 2 (([A , a y A z - a z A y ] + i{A x ,a y A y + a z A z }). (4) 

Here a = (6(t-t') sin[<j>(t) -cj)(t')}) p , C = (6(t-t')cos<j>(t'))p 
parametrize the evolution properties of the pulse at second 
order. The values of the parameters computed for some 
pulse shapes are listed in Tab. 1. 

Common pulse sequences. Transforming Eq. (4) appro- 
priately, we can now easily compute the result of application 
of any pulse sequence. In particular, the 7r-pulse X applied 
along the —x direction can be obtained from (—X) with a 
substitution a — > —a. As a result, e.g., the expansion of 
the evolution opeator for the one-dimensional sequence XX 
can be written as 

XX = 1 - 2tT p (Ao + a x A x ) - 2t 2 (A + o x A x f + C(t p 3 ), 

or it can be re-exponentiated as evolution with the effective 
Hamiltonian 

H xx = A ° + ^ + 0{r 2 p ). (5) 
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0.25 


Gooi 


0.0148978 


0.00735798 


0.249979 


G io 


0.148979 


0.0653938 


0.247905 


Si [4] 





0.0332661 


0.238227 


^2 [4] 





0.0250328 


0.241377 


Qi W 








0.239889 


Q2 [4] 








0.242205 



Table 1. Parameters of several common pulse shapes. The 
first line represents the "hard" 5-function pulse, Gooi denotes the 
Gaussian pulse with the width 0.01r p , while S n and Q„ denote 
the 1st and 2nd-order self-refocusing pulses from Ref. [4] . 



The corresponding calculation with the usual finite-width 
pulses (e.g., Gaussian) where s ^ produces a correction to 
the effective Hamiltonian already in the leading order, 

5H Tx = s{a z A y - a y A z ) + 0{s 2 t p ). 

This can be corrected by constructing a longer sequence, 
e.g., XXXX. Returning to pulses with s — 0, we list the 
expansions computed for several two-dimensional sequences: 

/ %GL 1 

Hxyxy = Ao + T p y — [A Z , A y ] - -[A ,a x A x - cr y A y ] 

~a v (Al + A 2 y ) + l±^a z {A x ,A y }) + 0(r 2 ), (6) 
h yxyxxyxy = ^o-^y (<Jy(A 2 x + Al) + i[A y ,A z ]j , (7) 

H YXYXXYXY = A + 0(t 2 ). (8) 

Atom in a cavity example. Consider an atom placed in a 
lossless cavity with a single resonant mode. The resonance 
part of the Hamiltonian can be written in the form (1), 
where A x = g(b + 6+), A y = ig(b - A z = 0, A = Atfb, 
g is a coupling constant, and A is the cavity frequency bias. 
Contrary to the conclusions of Ref. [7] , such a coupling can- 
not be suppressed with any one-dimensional pulse sequence. 

On the other hand, the two-dimensional four-pulse se- 
quence (6) provides a leading-order refocusing of the cou- 
pling. The subleading-order correction is present (the order 
of the sequence is K = 1), and it is not particularly small 
even for 2nd-order self-refocusing pulses with s = a = 0. 
The eight-pulse sequences (7) and (8) have equal accuracy 
with 2nd-order pulses but the 2nd-order accuracy of the lat- 
ter sequence is also retained with lst-order pulses. 

After tracing out the oscillator degrees of freedom, we 
can apply the results of Ref. [5] and expect the two 2nd- 
order sequences to provide an excellent refocusing accuracy 
even in an open system, as long as the refocusing rate is suf- 
ficiently high. We confirmed this expectation by a numerical 
simulation, where the bath was modeled as a classical cor- 
related random field. 

Conclusions. The main result of this work is the expan- 
sion (4) and the classification of the corresponding param- 
eters in Tab. 1. This allows an explicit computation of the 



error operators associated with refocusing in systems of ar- 
bitrary complexity. We illustrated the approach for several 
sequences applied to a qubit coupled to an oscillator. The 
8-pulse sequence (8) provides 2nd-order refocusing for any 
form of the coupling between the qubit and the oscillator. 
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